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Within an isospin-dependent transport model for nuclear reactions induced by 
neutron-rich nuclei, we perform a comparative study of ^^^Sn +^^^ Sn and 
^^^Zn+^^^Sn reactions at beam energies of 30 MeV/A and 400 MeV/A to identify 
optimal experimental conditions and observables for investigating the equation 
of state (EOS) of neutron-rich matter at RIA. Several observables known to be 
sensitive to the density- dependence of the symmetry energy are examined as 
a function of impact parameter for both reaction systems. In particular, the 
strength of isospin transport / diffusion is studied by using rapidity distributions 
of free nucleons and their isospin asymmetries. An approximate isospin equi- 
librium is established even for peripheral collisions with a very soft symmetry 
energy at 30 MeV/A. At 400 MeV/A, however, a strong isospin translucency is 
observed even in most central collisions with either soft or stiff symmetry en- 
ergies. Origins of these observations and their implications to determining the 
EOS of neutron- rich matter at the Rare Isotope Accelerator (RIA) are discussed. 
PACS numbers: 25.70.-z, 25.75.Ld., 24.10.Lx 
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The isospin-dependence of the nuclear equation of state (EOS) is one of the most im- 
portant but yet very poorly known properties of neutron- rich matter [1] . Its determination 
in laboratory- controlled experiments has profound implications to the study of many crit- 
ical issues in astrophysics [2,3]. Nuclear reactions induced by radioactive beams provide a 
unique opportunity to extract useful information about the EOS and novel properties of 
neutron-rich matter. In particular, at the Rare Isotope Accelerator (RIA), heavy-ions up to 
400 MeV/A will be available. It just enables the study of neutron- rich matter up to about 
3 times normal nuclear matter density. A number of experimental observables have already 
been identified as sensitive probes of the EOS of neutron-rich matter [4-17]. To gauge the 
sensitivity of some of these observables in reactions to be available at RIA, we carry out a 
comparative study of nuclear reactions nduced by ^^°Sn{N/Z — 1) and ^^^Zn{N/Z = 2.3) 
on ^"^^SniN/Z = 1.48) target at 30 MeV/A and 400 MeV/A from central to peripheral colli- 
sions. Results of this study are useful for planning experiments and designing new detectors 
at RIA. Moreover, the role of isospin degree of freedom in nuclear dynamics at RIA energies 
is an interesting subject in its own right. 

At present, nuclear many-body theories predict vastly different isospin dependence of the 
nuclear EOS depending on both the calculation techniques and the bare two-body and/or 
three-body interactions employed, see e.g., [18-21]. Various theoretical studies (e.g., [22,23]) 
have shown that the energy per nucleon e(p, S) in nuclear matter of density p and isospin 
asymmetry parameter 5 defined as 5 = iPn~ Pp)/iPn + Pp) can be approximated very well by 
a parabolic function e(p, 6) = e(p, 0) + Esym{p) ■ S"^. In the above e(p, 0) is the EOS of isospin 
symmetric nuclear matter and Esym{p) is the symmetry energy at density p. We shall use for 
isospin-symmetric nuclear matter a stiff EOS with Kq — 380 MeV which can reproduce the 
transverse flow data equally well as a momentum-dependent soft EOS with Kq = 210 MeV 
[24,25]. The choice of Kq does not affect our conclusions. The EOS of asymmetric matter 
should also be momentum-dependent since both the isoscalar and isovector potentials are 
momentum dependent [26,27]. For this comparative study, however, it is perfectly sufficient 
and much more efficient in calculations to use the momentum-independent EOS. 
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FIGURES 




FIG. 1. Symmetry energy (upper window) and symmetry potentials as a function of density 
for an isospin asymmetry of 5 = 0.2 and the 7 parameter of 0.5 and 2.0, respectively. 

The form of the symmetry energy as a function of density is rather strongly model 
dependent. We adopt here a parameterization used by Heiselberg and Hjorth- Jensen in 
their studies on neutron stars [28] Egymip) — Esym{po) ■ W^, where u = p/po is the reduced 
density and Esym{Po) is the symmetry energy at normal nuclear matter density po- By 
fitting the result of variational many-body calculations by Akmal et al [19], Heiselberg 
and Hjorth- Jensen found the values of EgymiPo) =32 MeV and 7 = 0.6. However, as 
shown by many other authors previously [20,29] using other approaches, the extracted value 
of 7 varies widely, even its sign at supranormal densities is undetermined. In a recent 
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analysis of isospin diffusion in fieavy-ion collisions at 50 MeV/A at the NSCL/MSU, some 
experimental indications of a 7 parameter as high as 2 was found [17]. Therefore, in this 
work we study influence of the symmetry energy on experimental observables by varying the 
7 parameter between 0.5 and 2. Assuming no momentum-dependence, the corresponding 
symmetry potentials can be derived directly from the density-dependent symmetry energy 
[11]. Shown in Fig. 1 are the symmetry energy (upper window) and symmetry potentials 
(lower window) at S—0.2 as a function of density. It is seen that the soft (7 = 0.5) symmetry 
energy leads to a larger (smaller) magnitude of the symmetry potential than the stiff one 
(7 = 2) at densities below (above) about po- Since the magnitude of the generally repulsive 
(attractive) symmetry potential for neutrons (protons) affects significantly the n/p ratio of 
emitted nucleons, the stiffness of the symmetry energy is thus expected to affect the isospin 
asymmetry of emitted free nucleons differently at low and high energies. To investigate the 
density dependence of the symmetry energy reactions with a broad range of beam energies 
are useful. 

Our study is based on an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) 
transport model (e.g., [11,4,30]). We select ^^^Sn and on the opposite boundaries 

of stabihty of mass number 100 as projectiles on a stable ^'^'^Sn target. The initial neutron 
and proton density distributions of the projectile nuclei were calculated by using the Hartree- 
Fock-Bogohubov method and were provided to us by J. Dobaczewski [31]. With these choices 
the ^^'^Sn -t--*^^^ Sn reaction is symmetric (highly asymmetric) in protons (neutrons) in the 
projectile and target, while the ^^^Zn 4-^^"^ Sn reaction is approximately symmetric (highly 
asymmetric) in neutrons (protons) . Comparing rapidity distributions of neutrons (protons) , 
for instance, of these two reactions will allow us to investigate the isospin transport and 
possible isospin equilibration. 

Shown in Fig. 2 and Fig. 3 are the distributions of free neutrons (left panels) and 
protons (right panels) as a function of scaled rapidity at a beam energy of 30 MeV/A and 
400 MeV/A, respectively, at an impact parameter of 1, 5 and 7 fm. Here we have selected 
free nucleons as those having local densities less than po/8 at freeze-out. 
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FIG. 2. Rapidity distributions of neutrons (left) and protons (right) in the reaction of 
1005.^ +124 and i°°Zn +^24 5-^^ at a beam energy of 30 MeV/A. 
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FIG. 3. Rapidity distributions of neutrons (left) and protons (right) in the reaction of 
1005.^ +124 sn and i°°Zn +^24 5-^^ at a beam energy of 400 MeV/A. 

First of all, it is seen that at 30 MeV/A the rapidity distributions of both neutrons and 
protons are about symmetric around {y /ybeam)cm = 0, indicating a high degree of isospin and 
kinetic equilibrium. At 400 MeV/A, however, distinct forward/backward asymmetries are 
seen in the rapidity distributions of neutrons in ^^^Sn +^^^ Sn and protons in ^^^Zn +^^^ Sn 
even in the most central collisions, indicating a significant nuclear translucency and non- 
stopping. Secondly, the main effect of the generally repulsive (attractive) symmetry potential 
for neutrons (protons) is to cause more (less) neutrons (protons) to become free. This 
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effect relies directly on the magnitude of the symmetry potential. At 30 MeV/A, since the 
symmetry potential is higher with 7 = 0.5 than with 7 = 2 as shown in Fig. 1, there are 
thus more (less) free neutrons (protons) emitted with 7 = 0.5 than with 7 = 2. While the 
opposite is true at 400 MeV/A since the magnitude of the symmetry potential reverses its 
order above po with 7 = 2 leading to higher symmetry potentials now. 




-0.2 I ■ ■ ' ■ ■ • ■ ■ ' ■ ■ • ■ ■ ' 1 

-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 



beam/cm 

FIG. 4. Rapidity distributions of isospin asymmetry of free nucleons in the reaction of 
Sn +^24 5^ Zn +^24 ^ heam energy of 30 MeV/A. 
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FIG. 5. Rapidity distributions of isospin asymmetry of free nucleons in the reaction of 
100 _^124 100 +124 ^ Y)eam energy of 400 MeV/A. 

To be more quantitative about the degree of isospin equilibrium and to examine effects 
of tfie symmetry potential on isospin transport, we show in Fig. 4 and Fig. 5 the isospin 
asymmetry of free nucleons as a function of rapidity. It is seen that a rather uniform isospin 
asymmetry is achieved over the whole rapidity range only with 7=0.5 at 30 MeV/A because 
of the large magnitude of the symmetry potential at densities relevant to the reaction. At 
400 MeV/A, an appreciable backward-forward asjTiimetric isospin asymmetry exists for both 
reactions. With 7 = 2 the overall isospin asymmetry of free nucleons are significantly higher 
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than that with 7 = 0.5. This is what one expects based on the density dependence of the 
symmetry potential shown in Fig. 1. Moreover, in contrast to the case at 30 MeV/A, the 
reaction time at 400 MeV/A is much too short compared to the time necessary for the 
system to reach the global isospin equilibrium [32,33]. 
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FIG. 6. Transverse momentum dependence of the neutron/proton ratio at mid-rapidity for the 
reaction of ^^^Sn +124 ^^^^ ^^^Zn +^^4 ^ beam energy of 30 MeV/A. 
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FIG. 7. Transverse momentum dependence of the neutron/proton ratio at mid-rapidity for the 
reaction of ^^^Sn -|-124 ^^^^ ^^^Zn +^^4 ^ beam energy of 400 MeV/A. 

A good quantitative measure of the isospin transport is the neutron to proton ratio at 
mid-rapidity. Shown in Fig. 6 and Fig. 7 are the n/p ratios at mid-rapidity versus the 
transverse momentum pt at 30 MeV/A and 400 MeV/A, respectively. For the neutron-rich 
system ^^^Zn +^^^ Sn at low transverse momenta at both beam energies, the n/p ratio is 
obviously increased compared to that of the reaction system. At both energies there is a 
general trend of an increasing n/p ratio at lower transverse momenta, except the case of 
7 = 0.5 at 30 MeV/A. This trend is mainly due to the Coulomb effect which shifts protons 
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from low energies to higher energies. While the generally attractive symmetry potential 
for protons works against the Coulomb potential. For protons at low transverse momenta 
the Coulomb potential dominates. While at higher transverse momenta with the parameter 
7 = 0.5 the symmetry potential dominates. With the soft symmetry energy at 30 MeV/A, 
the n/p ratio reaches the highest value at the highest transverse momentum. There is thus 
a strong sensitivity to the symmetry energy at high transverse momenta at 30 MeV/A. The 
n/p ratio differs by more than a factor of 2 by using 7 = 0.5 and 7 = 2. This is mainly 
due to the fact that these high transverse momentum nucleons have gone through regions 
of higher density gradients. 

At 400 MeV/A, effects of the symmetry potential is still obvious but not as dramatic as 
that at 30 MeV/A. Two effects are responsible for this observation. First, mean field effects 
are weaker at higher energies than at lower energies. Secondly, because of the particular 
functional form of the symmetry potential used here, the relative strengthes of the symmetry 
potentials are just opposite at supranormal and subnormal densities with 7 = 0.5 and 2 as 
shown in Fig. 1. At 30 MeV/A, only the low density part up to about 1.2po contributes 
while at 400 MeV/A a broad range of densities up to about 3po is involved. Of course, at 400 
MeV/A effects of the high density region dominate although there is a conceahng effect for 
particles going through both low and high density regions. Thus the n/p ratio is higher with 
7 = 2 than with 7 = 0.5 in contrast to the 30 MeV/A case. Therefore, the mid-rapidity n/p 
ratio observable is mostly useful for learning the symmetry energy at subnormal densities in 
reactions with beam energies around the Fermi energy. We shall discuss in the following two 
observables that are more useful for learning the symmetry energy at supranormal densities 
in reactions at higher beam energies. 



11 




FIG. 8. Neutron-proton differential flow for the reaction of Zn +1^^ Sn at a beam energy of 
30 MeV/A. 

We now examine the neutron-proton differential flow for the neutron-rich system 
^^^Zn +^^^ Sn. The neutron-proton differential flow in neutron-deficient reactions is not 
sensitive to the symmetry energy as one expects. The neutron-proton differential flow is 
defined as 

F:_^iy)^j:iptw,)/Niy), (1) 

i=l 

where Wi = 1(— 1) for neutrons (protons) and N{y) is the total number of free nucleons 
at rapidity y. It combines constructively effects of the symmetry potential on the isospin 
fractionation and the collective flow. It also has the advantage of maximizing effects of the 
symmetry potential while minimizing effects of the isoscalar potential [12]. Shown in Fig. 8 
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and Fig. 9 are the neutron-proton differential ffows in central and mid-central collisions at 
30 MeV/A and 400 MeV/A, respectively. 




FIG. 9. Neutron-proton differential flow for the reaction of Zn -j-^^** Sn at a beam energy of 
400 MeV/A. 

At 30 MeV/A, one observes an anti-flow while at 400 MeV/A there is a positive flow. 
This transition from negative to positive flow from the Fermi energy to relativistic energies 
is well known [34]. At both energies there is clearly an observable symmetry energy effects. 
Again, the higher magnitude of symmetry potential at low densities with 7 = 0.5 and at 
higher densities with 7 = 2 leads to the higher neutron-proton differential flow at 30 MeV/A 
and 400 MeV/A, respectively. Obviously the symmetry energy effect is now much stronger 
at 400 MeV/A than at 30 MeV/A. This is consistent with the well-known finding that the 
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transverse flow is sensitive to the early density and pressure gradients [3]. Therefore, the 
neutron-proton differential flow is more useful for studying the symmetry energy at high 
densities. 




FIG. 10. vr /vr^ ratio as a function of impact parameter in the reaction of ^^^Zn +1^^ Sn at a 
beam energy of 400 MeV/A. 

Another observable known to be sensitive to the high density behaviour of symmetry 
energy is the tt~/tt^ ratio [14,35]. Shown in Fig. 10 are the ti^/ti^ ratios at 400 MeV/A 
from central to peripheral collisions. It is seen that the ti~ /-k^ ratio is sensitive to both the 
n/p ratio of the reaction system and the symmetry energy, especially in central collisions. 
As discussed in detail in ref. [14,35], the 7r~/7r+ ratio is sensitive to the isospin asymmetry 
of the high density region {n/p) dense- While the latter is determined by the stiffness of the 
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symmetry energy. The soft symmetry energy at supranormal densities favors the formation 
of a neutron-rich high density region and a neutron-deficit low density region. In our dis- 
cussions above, it has been shown that at 400 MeV/A with 7 = 0.5 the n/p ratios of free 
nucleons are less than those with 7 = 2. Thus the corresponding {n/p)dense and tt'/tt'^ 
ratios are higher with 7 — 0.5. It is seen from Fig. 10 that the 7r~/7r+ ratio is about 25% 
higher with 7 = 0.5 than with 7 = 2 in central collisions. 

I. SUMMARY 

In summary, within an isospin-dependent transport model for nuclear reactions induced 
by neutron-rich nuclei, we carried out a comparative study of ^'^^Sn+^'^^Sn and ^^'^Zn+^^'^Sn 
reactions at a beam energy of 30 MeV/A and 400 MeV/A, respectively. The strength of 
isospin transport /diffusion is studied by using the rapidity distributions of free nucleons 
and their isospin asymmetries. An approximate isospin equilibrium can be obtained even 
for peripheral coUisions with only the soft symmetry energy at 30 MeV/A. At 400 MeV/A, 
however, an appreciable isospin translucency happens even in the most central collisions. 
We found that the neutron/proton ratio at mid-rapidity but high transverse momenta at 30 
MeV/A is most sensitive to the symmetry energy at subnormal densities. On the other hand, 
at higher beam energies, e.g., 400 MeV/A, neutron-proton differential flow and the 7r~/7r+ 
ratio are found to be more sensitive to the density- dependence of the symmetry energy at 
supranormal densities. To map out the entire density-dependence of the symmetry energy 
a combination of several complementary observables in a broad range of beam energies 
is necessary. These observations are useful for planning experiments and designing new 
detectors for exploring the EOS of neutron-rich matter at RIA. Several observables useful 
for studying the EOS of neutron-rich matter require the detection of neutrons. 

We would like to thank Dr. J. Dobaczewski and Dr. Lie- Wen Chen for providing us the 
initial nucleon density distributions used in this work. This work was supported in part by 
the National Science Foundation Grant No. PHY-0088934 and PHY-0243571. 
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